Interplay between valence and core excitation mechanisms in the breakup of halo 

nuclei 
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The phenomenon of core excitation in the breakup of a two-body halo nucleus is investigated. We 
show that this effect plays a significant role in the reaction dynamics and, furthermore, its inter- 
ference with the valence excitation mechanism has sizable and measurable effects on the breakup 
angular distributions. These effects have been studied in the resonant breakup of 11 Be on a carbon 
target, populating the resonances at 1.78 MeV (5/2+) and 3.41 MeV (3/2 + ). The calculations have 
been performed using a recently extension of the DWBA method, which takes into account the effect 
of core excitation in both the structure of the halo nucleus and in the reaction mechanism. The 
calculated angular distributions have been compared with the available data [Fukuda et al, Phys. 
Rev. C70, 054606]. Although each of these resonances is dominated by one of the two considered 
mechanisms, the angular patterns of these resonances depend in a very delicate way on the inter- 
ference between them. This is the first clear evidence of this effect but the phenomenon is likely to 
occur in other similar reactions. 

PACS numbers: 24.50.+g, 25.60.Gc,27.20.+n, 25.40.Ep 
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Introduction - The study of exotic nuclei has played a 
key role in nuclear physics over the past 25 years. Our 
current knowledge of their peculiar structural proper- 
ties comes mainly from measurements of removal cross 
sections, transfer and breakup reactions. Breakup re- 
actions have provided useful information on the ground 
state properties, such as binding energies, spectroscopic 
factors and angular momentum (eg. [l|, [2j]). Moreover, 
when exclusive measurements are possible, i.e., all out- 
going fragments are detected after breakup, these experi- 
ments can be used to infer spectroscopic properties of the 
continuum, such as the location and spin/assignment of 
resonant states and dipole strengths 0, 00|- 

In the case of halo nuclei, loosely bound exotic nuclei 
composed of a tightly bound core surrounded by one or 
two loosely bound nucleons, these processes have been 
conveniently modeled using a three-body model, com- 
prising the two-body weakly bound projectile and the 
target. This has motivated the development and revival 
of few-body theories, such as the Continuum-Discretized 
Coupled-Channels (CDCC) method the adiabatic 
(frozen-halo) approximation @-ll|, a variety of semiclas- 
sical approaches fl2T4l5| and, more recently, the Faddeev 
equations 

In their standard formulations, these methods assume 
a single-particle description of the valence particle rela- 
tive to the core. Possible excitations of the core are ne- 
glected or, at most, taken into account effectively through 
the core-target optical potential. Although this approach 
has been used with relative success in the analysis of 
many reactions, it has been recently shown 2(| 21 1 that 
this simplified picture is not always accurate, due to the 
effects of core excitation. This affects in two ways. First, 
the presence of core admixtures in the states of the pro- 
jectile means that these states cannot be simply treated 



as single-particle states calculated in some mean field po- 
tential. Second, the interaction of the core with the tar- 
get may give rise to transitions between these core states, 
leading also to the breakup of the projectile. Although 
these two effects have been commonly ignored in the 
analysis of reactions with exotic beams, some progress 
has been made in recent years toward their incorpora- 
tion in existing reaction formalisms. For example, Sum- 
mers et al. [22| have proposed an extended version of 
the CDCC method which treats the structure of the few- 
body projectile within the particle-rotor model. More 
recently, a simple extension of the standard DWBA am- 
plitude which takes into account these effects approxi- 
mately has been proposed [2(| 21 1 . Using this method, it 



was found that dynamic core excitation are indeed very 
important to describe the breakup of n Be on protons at 
^70 MeV/nucleon. Some effects have been also found in 
the elastic scattering of 8 B on a carbon target |23|. using 
an extension of the adiabatic model of Ref. Ilfj. 
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Although the calculations presented in [2C 
vide a clear evidence of the importance of core excita- 
tion in nuclear breakup, the restricted energy and angu- 
lar resolution of the analysed data prevented a detailed 
assessment of the relative importance and the interplay 
between the valence and core excitation mechanisms. In 
this Letter, we present new results showing that the pres- 
ence of core admixtures in the halo nucleus and the sub- 
sequent dynamic core transitions give rise to very distinc- 
tive effects on the shape and magnitude of the breakup 
cross sections angular distributions. Moreover, it is found 
that these effects depend very critically on the amount of 
core excitation of the halo nucleus. This sensitivity opens 
new possibilities for extracting spectroscopic information 
of halo and other weakly bound nuclei, by comparing the 
measured angular distributions with a reliable reaction 



2 



model. To illustrate these effects, we present here calcu- 
lations for the breakup of the one-neutron halo nucleus 
11 Be on a carbon target at a bombarding energy of 70 
MeV/nucleon in particular to the excitation of the 
resonances at E x = 1.78 MeV and 3.41 MeV. The calcu- 
lations are performed using the extended version of the 
DWBA method of Refs. El 
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Reaction model- We now describe briefly the core- 
excitation model used here. We outline here the main 
formulas, and refer the reader to Refs. [20I |2~0 ] for fur- 
ther details. We consider the inelastic excitation of a 

to 
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projectile nucleus, initially in its ground state, ^ 
a state ^ji M i (bound or unbound). Within a two-body 
(core+valence) description of the projectile, these wave- 
functions are expanded as 



JM 



(1) 



where the functions <p a (r) describe the relative motion 
between the valence particle and the core and $jm c (£) 
are the core eigenstates with angular momentum / and 
projection M c . The index a denotes the set of quantum 
numbers {£, j, I}, with £, s and j being the orbital angular 
momentum, the intrinsic spin of the valence particle and 
their sum (j = £ + s), respectively. The functions (p a (r) 
and <&iM c (0 depend on the assumed structure model, to 
be specified later. 

Consistently with the assumed two-body description of 
the projectile, the transition potential for this breakup 
process is the sum of the valence-target and core-target 
interactions, i.e., Vt = V vt (R v t) + V ct (R c t,£,)- While 
the valence-target interaction is taken to be central and 
to depend exclusively on the valence-target separation, 
the core-target interaction is assumed to depend on the 
core internal degrees of freedom, and can therefore induce 
transitions between different core states. By expanding 
this interaction in multipoles (A) and separating the cen- 
tral (A=0) part, the DWBA amplitude, describing the 
excitation of the halo nucleus during the collision, splits 
into two terms: 
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where K (K 1 ) is the initial (final) linear momentum and 

x\ +) {R) (X { f + \R)) the initial (final) distorted wave de- 
scribing the projectile-target relative motion. The transi- 
tion amplitude given by Eq. (j3]) (valence amplitude here- 
after) describes excitations between different valence con- 
figurations, but without altering the state of the core. 
This term is evaluated following the standard techniques 



used in coupled-channels codes. The second term in 
([5]) , denoted core excitation amplitude, contains the non- 
central part (A>0) of the core-target interaction and can 
therefore produce core transitions. Consequently, this 
term accounts for the dynamic excitation of the core dur- 
ing the collision. In [20i . |2~0 ] , it was shown that this term 
acquires a very simple form when evaluated in the no- 
recoil approximation (i? ct fts R), 
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is a geometric factor [20j, [21 1, R^ are 



the radial parts of the <p a (r*) functions and T C { is re- 
lated to the core-target two-body transition amplitude 
for a core transition IM C — > I'M' C of multipolarity A as 



T, 
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= (IM c \Li\rM' c )T^ } 
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Results. - The model has been applied to the resonant 
breakup of n Be on a 12 C target at 70 MeV/nucleon, mea- 
sured at the RIKEN facility by Fukuda el al. [4]. The 
relative energy spectrum of the 10 Be+n center of mass 
shows peaks at E x =1.78 MeV and £' x =3.41 MeV. Their 
angular distributions were compared with DWBA calcu- 
lations, based on the vibrational collective model, sug- 
gesting a A = 2 transition for both states. These states 
were identified with the 5/2^ and 3/2^ resonances pre- 
dicted by shell-model calculations. 

In the present work, we reanalyze these data using 
the aforementioned core-excitation DWBA model. The 
structure of the 11 Be nucleus is described within the 
particle-rotor model (PRM) of Bohr and Mottelson, with 
the parameters given by the model Bel lb of Ref. 
This model assumes a permanent quadrupole deforma- 
tion for the 10 Be core with ,02=0. 67. 

The 11 Be wave functions are obtained by diagonal- 
izing the 11 Be Hamiltonian in a particle+core basis of 
the form !$/(£) ® V™? M (f))jM I where ^™°,(r) (with 
n=l,. . . ,N) are a truncated set of Transformed Harmonic 
Oscillator (THO) functions [2^], which are used as a basis 
for the valence-target relative motion. This basis is ob- 
tained by applying an analytic local scale transformation 
(LST) to the conventional harmonic oscillator basis. The 
modelspace is restricted to I 7 ' — + ,2 + and £ < 2. For 
the ground state, a basis of N=15 oscillator functions was 
used. Resonant states are identified with stabilized en- 
ergies with respect to the basis size (N). The parameters 
of the LST are the same as those used in [25(. 

The components involved in our calculations allowed 
by the present model space, along with their respective 
weights (spectroscopic factors) are listed in Table H Also 
listed in this table are the shell-model spectroscopic fac- 
tors obtained with the code OXBASH, using the WBT 
interaction of Warburton and Brown [26] . Both models 
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predict very similar spectroscopic factors for the ground 
state and the 5/2+ resonance, and only some small dif- 
ferences are found in the 3/2 + state. The ground state 
corresponds predominantly to a | 10 Be(0 + ) ® S1/2) config- 
uration, with some admixture of the | 10 Be(2+) ® ^5/2) 
configuration. The 5/2 + state is mainly based on the 
10 Be ground state. On the other hand, the 3/2+ reso- 
nance is mainly built on top of the excited core. Accord- 
ing to this result, it is expected that the population of 
the 5/2 + state is mainly due to the valence excitation 
mechanism, whereas the excitation of the 3/2+ state will 
be mostly due to a core-excitation mechanism. 

To illustrate the sensitivity of the calculation with the 
structure model, we have considered two additional mod- 
els assuming pure single-particle configurations for the 
n Be g.s. and the 5/2 + and 3/2+ resonances. For the 
11 Be(g.s.) we consider a pure |0 + ®2s 1 / 2 ) configuration. 
For the 5/2 + resonance we consider two single-particle 
models: i) |0+ ® ld 5 / 2 ) (denoted SP1) and ii) \2 + ®2s 1/2 ) 
(SP2). In the former, the resonance is populated by 
means of a valence excitation mechanism, whereas in the 
second model the excitation is due to a pure core excita- 
tion effect. Similarly, for the 3/2 + we consider also two 
extreme models: i) |0 + ®ld 3 / 2 ) (SP1) and ii) 2+®2s 1/2 ) 
(SP2). The required radial wavefunctions are taken from 
the PRM calculation, conveniently normalized to one. 

The n+ 12 C potential was taken from Ref. |27j. The 
central and transition components of the 10 Be+ 12 C po- 
tential were generated by a double folding procedure, 
convoluting an effective nucleon-nucleon (NN) interac- 
tion with the 10 Be and 12 C matter densities. The latter 
were taken, respectively, from the Antisymetrized Molec- 
ular Dynamics (AMD) calculation of Ref. [28| and from 
the parametrization of Ref. [29] . For the effective NN in- 
teraction we adopt the spin-isospin independent part of 
the M3Y interaction (3^| based on the Reid soft-core NN 
potential. For the imaginary part of the 10 Be+ 12 C po- 
tential we assume the same geometry as for the real part. 
A renormalization factor was included to reproduce the 
elastic scattering data of 10 Be+ 12 C at 59.4 MeV/nucleon 
from Ref. [3l|. Further details of these calculations will 
be provided elsewhere. 

In Fig. Q] we compare the calculated angular distribu- 
tions with the experimental data of Ref. [4]. The up- 
per and bottom panels correspond to the 5/2+ (E x =1.78 
MeV) and 3/2+ (E X =3A1 MeV) resonances. It is readily 
seen that the pure single-particle models SP1 and SP2 do 
not reproduce the shape of the resonances. In the model 
SP1 (pure valence excitation) the maxima and minima 
are shifted to smaller angles with respect to the data and 
the angular distribution decay too fast. On the other 
hand, in the model SP2 (pure core excitation mechanism) 
the maxima and minima are shifted to larger angles. Fi- 
nally, the full PRM model, which includes both valence 
and core excitation mechanisms and their interference, 
the position of the maxima and minima is very well re- 
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FIG. 1. (Color online) Angular distribution for the E x = 1.78 
MeV and 3.41 MeV states in lx Be. The circles are the data 
from Ref. 0]. The curves correspond to the extended DWBA 
calculations, including core excitation effects, using different 
structure models for the 1J Be nucleus. For the single-particle 
models (SP1 and SP2) the resonance configuration is indi- 
cated in the labels. 



produced. It is also seen that the absolute magnitude of 
the data is overestimated. Except for this discrepancy 
in the normalization, it is clear that the shape is appre- 
ciably improved with respect to the pure single-particle 
description and that the interference between the valence 
and core excitation mechanisms is crucial to account for 
the correct shape of the oscillations. 

It is enlightening to consider separately the contri- 
bution of the valence and core excitation amplitudes, 
Eqs. © and (gj). These are depicted in Fig. [5] for the 
PRM model. In this plot, the calculations have been 
convoluted with the experimental angular resolution Q 
for a more meaningful comparison with the data. As 
anticipated, the 5/2 + resonance is mainly populated by 
the valence excitation mechanism, due to its dominant 
lo Be(0 + ) configuration, whereas for the 3/2 + state the 
dynamic core excitation mechanism is the dominant one. 
It is also seen that both contributions are out of phase, 
and the interference between them is very important. In 
fact, none of them separately is able to reproduce by it- 
self the position of the maxima and minima of the data, 
whereas their coherent sum (solid line) reproduces very 
well this pattern. This result illustrates very nicely the 
delicate interplay between the valence and core excitation 
mechanisms in the breakup of a deformed halo nucleus, 



4 



TABLE I. Spectroscopic factors for the ground state and resonant wavefunctions of 11 Be, according to the particle-rotor model 
(PRM) and the shell-model calculations (WBT) presented in this work. 





State 


Model 


|0+ ® (£s)j) 


|2+®s 1/2 ) 


|2+®d 3/2 ) 


|2+®d 5/2 ) 




1/2+ (g.s.) 


PRM 


0.857 




0.021 


0.121 






WBT 


0.762 




0.002 


0.184 


5/2+ 


(E x =1.78 MeV) 


PRM 


0.702 


0.177 


0.009 


0.112 






WBT 


0.682 


0.177 


0.009 


0.095 


3/2+ 


(£^=3.41 MeV) 


PRM 


0.165 


0.737 


0.017 


0.081 






WBT 


0.068 


0.534 


0.008 


0.167 
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FIG. 2. (Color online) Valence (dashed) and core (dot- 
dashed) contributions to the breakup of the 1.78 and 3.41 
MeV resonances populated in the 11 Be+ 12 C reaction at 70 
MeV/nucleon, using a particle-core description of the n Be 
nucleus. The solid is the coherent sum of both contributions. 



in the breakup of halo nuclei using and recently proposed 
extension of the DWBA method. We have shown that the 
presence of core admixtures in the initial and final states 
has a sizable impact in the interference pattern of the 
breakup cross section and hence a high sensitivity on the 
underlying structure model of the halo nucleus. This ef- 
fect has been evidenced for the first time in the scattering 
of n Be on 12 C at 70 MeV/nucleon, where we have shown 
that the inclusion of these core excitation effects improves 
significantly the agreement with the data 0| and provides 
very valuable spectroscopic information, which would be 
very difficult to extract from other methods. Finally, we 
emphasize that, although the calculations have been pre- 
sented for the u Be nucleus, we do expect these effects 
to be important in other relevant cases, such as in the 
breakup of the odd carbon isotopes 15 < 17 < 19 C. 
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ing us the 10 Be microscopic densities and to T. Naka- 
mura for his help regarding the 11 Be+ 12 C data and 
the convolution with the experimental resolution. This 
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isterio de Ciencia e Innovation under project FPA2009- 
07653, and by the Spanish Consolider-Ingenio 2010 Pro- 
gramme CPAN (CSD2007-00042). J.A.L. acknowledges a 
research grant by the Ministerio de Ciencia e Innovation. 



like 11 Be. Note that the weak contribution of the valence 
mechanism in the 3/2 + case is a consequence of the small 
spectroscopic factor associated to the |0 + ®d 3 / 2 ) config- 
uration (see Table fl} . This fact explains also that this 
resonance is very weakly populated in transfer reactions, 
such as 10 Be(c?,p) 11 Be [32j, making it difficult the extrac- 
tion of spectroscopic information from these experiments. 
In these cases, the approach presented in this work, based 
on the analysis of breakup reactions, provides a powerful 
alternative to access this information. 

Conclusions - In conclusion, we have studied the inter- 
play between the valence and core excitation mechanisms 
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